Phonon Hall Effect in Nanoscale Four-Terminal Junctions 
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Using an exact nonequilibrium Green's function formulism, the phonon Hall effect (PHE) for 
paramagnetic dielectrics is studied in a nanoscale four-terminal device setting. The temperature 
difference in the transverse direction of the heat current is calculated for two-dimensional models 
with the magnetic field perpendicular to the plane. We find that there is a PHE in nanoscale 
paramagnetic dielectrics, the magnitude of which is comparable to the millimeter scale experiments. 
If the dynamic matrix of the system satisfies mirror refiection symmetry, the PHE disappears. The 
Hall temperature difference changes sign if the magnetic field is sufficiently large or if the size 
increases. 

PACS numbers: 66.70.-f, 72.15.Gd, 72.20.Pa 
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In parallel to the study of electronics, phononics gets 
very active and hot recently Several conceptual 
phononic (thermal) devices such as thermal diode, tran- 
sistor, logic gate and even thermal memory have been 
proposed to control phonon and process information with 
phonons. The magnetic field is another degree of freedom 
could be potentially used to control phonon transport Q • 
Indeed, a novel phenomenon ~ phonon Hall effect (PHE) 
- has been discovered experimentally by Strohm, et al. 
Q , which is an analog of the electrical Hall effect for the 
heat flow in dielectrics. The authors found a tempera- 
ture difference up to 200 /iK between the sample edges 
in the direction perpendicular to both the heat flow and 
the magnetic field. This effect has been confirmed in 
The electronic Hall effect is well understood in terms of 
Lorentz force. However, since no charges are involved and 
phonons can not couple to the magnetic field directly, it 
is thus not trial and straight forward to understand the 
PHE. The magnetic field can polarize the paramagnetic 
ions; the subsystem of isolated ions carrying magnetic 
moment M couples to phonons, this spin-phonon inter- 
action (SPI) determines the phonon Hall effect [1, 01- 

Theoretical models for PHE have been proposed in 
Refs. [1, 0|, in which phonons are treated ballistically. 
However, according to Ref. jj], the mean free path (1/im) 
is far less than the system size (15.7 mm), therefore, it 
is not appropriate to treat the diffusive PHE with bal- 
listic theory. Moreover, in all the previous theoretical 
work, the spin-phonon interactions were considered by 
perturbation theory, and they are not consistent with 
each other. 

In this Letter, we would like to address two questions: 
(1) whether the PHE exists in nano-scale systems, or in 
other words whether the diffusive PHE can be observed 
in the ballistic regime; (2)treat the PHE by an exact - 
non-perturbative theory. The first question is in fact not 
trivial at all. Since many physics laws valid in macro- 
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FIG. 1: (color online) The four-terminal PHE setup used for 
calculating the thermal conductance and the temperature dif- 
ference T3 — T4,. (a) The left and right leads have temperatures 
Ti and T2, the upper and lower probe- leads have temperatures 
T3 and Ti. The center part can be different lattices, such as 
square lattice (b) or honeycomb lattice (c). 



scopic scale are not necessary true in nanoscale. For ex- 
ample, Fourier law of heat conduction is broken down in 
nanoscaleQ- As for the second question, we will take 
the nonequilibrium Green's function (NEGF) approach. 
The NEGF is an elegant and powerful method to treat 
nonequilibrium and interacting systems in a rigorous way. 
NEGF is widely applied to the electronic and thermal 
transport, and is successful to study the spin Hall effect 
in junctions [8]. 

To develop a nonperturbative theory for PHE in 
nanoscale four terminal junctions, we consider a model 
as shown Fig. [T] by taking into account the actual mea- 
suring process. The thermal conductance of the system 
can be calculated by the NEGF method. As we shall see 
later, our model systems can produce features similar to 
experiments, even though our systems are of nanometer 
scale while the experimental systems are of millimeter 
scale and are in the diffusive regime. 

We consider the same Hamiltonian of SPI as in 
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Refs. 0, @], which can be expressed in the form Hi = 
gJ2^n ■ {Un X Pn). Here, Un and Pn are the vectors of 

n 

displacement and momentum of the n-th lattice site. In 
the presence of a magnetic field B, each lattice site has 
a magnetization M. For isotropic SPI, the isospin s„ is 
parallel to M„, and the ensemble average of the isospin 
is proportional to the magnetization, that is (sVi) = cA/. 
Therefore, under the mean-field approximation, the SPI 
can be represented as 



Hi=J2^-iUnXP.,,), 



(1) 



where, A = gcM, has the units of frequency. According 
to 0, A is estimated to be 0.1 cm^^ w 3 x 10^ Hz at 
B = IT and T = 5.45 K, which is within the possible 
range of the coupling strength in ionic insulators [lo| . 
In our calculation, we will use this relation to map A to 
magnetic field. 

Now we turn to the derivation of the general formula 
for the temperature difference in the system as illustrated 
in Fig.[l] where a 2D lattice sample, which can be honey- 
comb lattice or square lattice, is connected with four ideal 
semi- infinite leads. We denote the lattice as Nfi x Nc, 
A^-Ri correspond to the number of rows and columns. 
We adjust the temperatures of upper and lower probes 
and T/^ such that the heat currents from these two leads 
vanish, namely, = I/^ = 0. Then we get the relative 
Hall temperature difference as i? = (Tg - Ti)/{T^ T2). 

We treat the PHE problem by NEGF method [l3] anal- 
ogous to those used in thermal transport [l4'|. The total 
Hamiltonian is assumed to be 
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U^APo, 



where Ha = \ [PaPa + U^K^Ua), and A is an antisym- 
metric, block diagonal matrix with the diagonal elements 

^. Here, the integers to 4 are associated with 

-A y ' ^ 

the center region, left, right, upper, and lower leads, re- 
spectively. Ua (Pa) are column vectors consisting of all 
the displacement (momentum) variables in region a. Ka 
is the spring constant matrix and Vpfi = (Vo./s)^ is the 
coupling matrix between the /3 lead and the central re- 
gion. The dynamic matrix of the full linear system with- 
out SPI is 



K = 



( 
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(3) 



We obtain the equation for [/o and Pq as 



dr 



Po{t) ~ AUo{t), (4) 



= -KoUoir) ~ J2 VoMt) - APo{t). (5) 

/3=1 

The energy flux to the central region from the lead a is. 



Hr, 



{[Ha,H]), a = 1,2,3,4. (6) 



We define the contour-ordered Green's function as 
G"'^(r,r') = -{ {T^UaiT)Up{T'f) , where a and (3 re- 
fer to the region that the coordinates belong to and 
Tc is the contour-ordering operator. Then the equa- 
tions of motion of the contour ordered Green's func- 
tion can be derived. In particular, the retarded Green's 
function for the central region in frequency domain is 



G''M= {uj + irj)'^ - Ko-J:''[uj]- + 2iujA 



Here, 



S*" = ^ EJ^, and Sq, = Vo,q5qKi,o is the self-energy due 

to interaction with the heat bath, = [{uj+irj)'^ — Ka]~'^ ■ 
The lesser Green's function is obtained through — 
QrY,<Q"- in the usual way. We thus can calculate the 
heat flux by the following formula. 
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(7) 



If the temperature differences among the leads are very 
small, we can regard the system as at linear response 
regime, — T + A^- The linearized heat flux from each 
heat bath can be written as 



(8) 



13=1 



(2) The conductance from heat bath a to ^ is deflned as 



(9) 



where, T^j^M = TriG^TpG^Ta), f = (e^-/^^^ - l) \ 
and Fq = z(S^[u;] — [a;]). Therefore, if we set the heat 
flux of lead 3 and lead 4 to zero, we can get the relative 
Hall temperature difference as 



R = (A3-A4)/(Ai-A2) 

0'130'24 — cr23Cri4 



(cri3 + (T23 + 0'43)(o-i4 -t- (724 + CT34) — a^^a^^ 



.(10) 



Equations ([8|) and ([9]) are the Landauer-Biittiker theory 
11, l2| applied to the multiple-lead thermal transport. 

In the following calculation. We assume a lattice 
constant a = 2.465 A, and the force constant — 
0.02394 eV/(amu-A2), Kt = KL/'i- The ratio of the lon- 
gitudinal and transverse sound speed to he S = vl/vt 
^/KlJKt ~ 2. Then the speed of sound for longi- 
tudinal acoustic phonons is about 4000 m/s. As men- 
tioned above, A is estimated to be about 3 x 10^ Hz « 
2.0 X 10-6 eV at B = IT. 
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FIG. 2; (color online) Hall temperature difference R versus 
magnetic field B at temperature T = 5.45 K. The hexagon 
and square line correspond central regions for the honeycomb 
and square lattices with a nearest-neighbor coupling. The red 
dotted line corresponds to a linear fit from to 40 T. The size 
of the center region for honeycomb lattice is 9 x 6, the same 
with the inset (c) in Fig. [1] 



FIG. 3: (color online) Thermal conductance versus the mag- 
netic field at temperature T = 5.45K for the honeycomb lat- 
tice, (a) shows the conductance between two longitudinal 
leads (Ti2. (b) shows the conductance between one longitudi- 
nal lead and one transverse probe-lead. The circle and trian- 
gular lines correspond to ais and ai4, respectively. The size 
of center region is 9 x 6. 



From Eq. (|10p . we find that the relative HaU temper- 
ature difference R is an odd function of magnetic field. 
Onsager relation, cra^(A) — cr^c[(— A), always holds due 
to the definition of the conductances. Furthermore, if 
there is a symmetry operation S such that 

SKS-^^K, SAS-'^ = -A, (11) 

then, (Tq^(A) — crQ^(— A). If this relation is true, then 
there is no phonon Hall effect in the system. These sym- 
metry considerations are consistent with a different treat- 
ment for bulk systems based on Green-Kubo formula [3] ■ 

We discuss numerical results in the following. Fig. [2] 
shows the temperature difference changing with magnetic 
field at temperature T = 5.45 K for the honeycomb and 
square lattices with nearest-neighbor couplings. For the 
honeycomb case, the Hall temperature is odd and linear 
in the magnetic field between and 40 T, in that range 
the slope of the curve is 3 x 10~^K/T, comparable to 
the experimental data in Rcf. When the magnetic 
field is extremely large, it will decrease. From our cal- 
culation, we find that the triangular lattice has a similar 
behavior. However, for square lattice with the nearest- 
neighbor coupling, there is no PHE at all. The spring 
constant matrix between every nearest coupling sites is 
diagonal for the square lattice. This matrix and also the 
full matrix K is invariant with respect to a refiection in 
X OT y direction, thus satisfying Eq. (jlip . If we consider 
next-neighbor couplings of the lattice, the dynamic ma- 
trix K will not have the mirror reflection symmetry, and 
the PHE can come out. 

We show the conductances among different leads in 
Fig. [31 We set all the couplings between the leads and 
central region the same, and all the leads and central 
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FIG. 4: (color online) The relative Hall temperature difference 
R versus the large magnetic field B (a) and high equilibrium 
temperature T for honeycomb lattice (b). (a) square shows R 
changing with the magnetic field (left scale) , the red solid line 
shows the conductance difference (J13 — ai4 versus magnetic 
field (right scale), (b), R vs. equilibrium temperature at 
B = IT. 



region have the same spring constants for simplicity. Be- 
cause of the symmetry of the system, we have additional 
relations, ctis = 0-32 = (T24 = 0-41, and an = (742 = 0-23 = 
(T31. We find that the conductance between two longi- 
tudinal leads or two transverse probe-leads are even in 
the magnetic field, which can be seen in Fig. [3{a), (734 
has the same property. However, for honeycomb lattice 
the conductance between one longitudinal lead and one 
transverse probe-lead is not an even function of mag- 
netic field [Fig. mb)], which gives contribution to the 
Hall temperature difference. Therefore, for honeycomb 
lattice, the temperature difference is not zero. But for 
square lattice, (T13 is an even function of magnetic field, 
the same is true for other components, no PHE exists in 
such systems. 
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FIG. 5: (color online) (a) The relative Hall temperature dif- 
ference versus magnetic field for difi'erent ratio of the longi- 
tudinal and transverse sound speed 5 = vl/vt- The inset 
shows R versus J at B = IT. The data are calculated for 
9x6 honeycomb lattice at T = 5.45K. (b) The relative Hall 
temperature difference versus the number of rows of atoms for 
fixed aspect ratio Nr : Nc = I : 2 at B=l T and T = 5.45 
K. 



We show the numerical results for the ratio R at 
T = 5.45K for honeycomb lattice in Fig. IDJa). The tem- 
perature difference will not be linear when the magnetic 
field is larger than 40 T, after about HOT, it will de- 
crease, and about Be ~ 230 T, R changes sign to neg- 
ative. It is the same critical point for the difference of 
conductances cria — (T14, which is consistent with Eq. (jlOp . 
In Fig. m^b), we show R versus temperature at B = 1 T. 
When the temperature increases, R will increase almost 
linearly. After some value, it decreases, and then in- 
creases again. At last, it tends to a constant. This be- 
havior is due to the competition of the numerator and de- 
nominator in Eq. (llOp . When temperature is very high, 
all the conductances tend to constants due to ballistic 
thermal transport. 

In Ref . [5] , it is reported that R decreases with increas- 
ing ratio of the longitudinal and transverse sound speed 
S = vl/vt and changes sign when 5 becomes large than 
5. However, we find that when the ratio {5 > 1) become 
large, R increases, see Fig.[5l[a). At exactly (5=1, when 
the longitudinal speed equals to the transverse speed, 
there is no PHE, which testifies our condition, Eq. pT|) . 
for the absence of PHE. All the spring constant matrices 
between the nearest- neighbors become diagonal &i 5 — 1, 
the condition Eq. pip holds for a mirror reflection op- 
eration. If (5 < 1, i? increases again with the decreasing 
of 5. Although the ratio R does not change sign with (5, 
due to the ballistic nature of a small system, the ratio 
R is sensitive to the geometric details, which is shown in 
Fig. mb) , the magnitude and the sign of R will change 
with the size increasing. 

In conclusion, a theory for PHE in nanoscale param- 
agnetic dielectrics by NEGF approach is developed. The 
results are consistent with the essential experimental fea- 
tures of PHE, such as the magnitude and linear magnetic 
field dependence of the observed transverse temperature 



difference. We find that there is no PHE if the lattice 
satisfies a certain symmetry. The symmetry of the dy- 
namic matrix K is the key point for the existence of 
PHE. The Hall temperature difference changes with equi- 
librium temperature and tends to be a constant at last. 
And the Hall temperature difference does not change sign 
with the ratio of the longitudinal and transverse sound 
speed in the range of (5 € (0.1, 10). 

Most of our results should be verified by experiments 
on nano-scale paramagnetic dielectrics, which can have 
potential applications to controlling nanoscale phonon 
transport. For most paramagnetic dielectric materials, 
because of the complexity of the coupling (such as next 
or next-next neighboring interaction), the dynamic ma- 
trix does not satisfy the mirror refiection symmetry, the 
PHE can be present. The Hall temperature difference be- 
havior with magnetic field can be measured in a strong 
magnetic field (if A = 2.0 x 10~^ eV, the magnetic field 
will be ten times smaller). From our study, the PHE can 
be measured in nanoscale system at relative high tem- 
perature 100 K). 
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